Scaling of cluster heterogeneity in percolation transitions 
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We investigate a critical scaling law for the cluster heterogeneity H in site and bond percolations 
in d-dimensional lattices with d = 2, • • • , 6. The cluster heterogeneity is defined as the number of 
distinct cluster sizes. As an occupation probability p increases, the cluster size distribution evolves 
from a monodisperse distribution to a polydisperse one in the subcritical phase, and back to a 
monodisperse one in the supercritical phase. We show analytically that H diverges algebraically 
approaching the percolation critical point p c as H ~ \p— p c |~ with the critical exponent a associ- 
ated with the characteristic cluster size. Interestingly, its finite-size-scaling behavior is governed by 
a new exponent vh = /d)v where df is the fractal dimension of the critical percolating cluster 

and v is the correlation length exponent. The corresponding scaling variable defines a singular path 
to the critical point. All results are confirmed by numerical simulations. 
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I. INTRODUCTION 

Percolation is a geometric phase transition for connec- 
tivity [l|, 0] . Suppose that a fraction p of sites or bonds 
are occupied in an infinite lattice. When p is less than a 
threshold p c , any site is connected to others up to a fi- 
nite distance via occupied sites or bonds. As p increases, 
clusters (sets of connected sites) grow until there emerges 
a spanning, giant, or an infinite cluster to which a finite 
fraction of sites belong. 

The percolation transition is manifested in several 
quantities. The order parameter to, given by the fraction 
of sites belonging to the giant cluster, becomes nonzero 
beyond the transition. The mean cluster size S, defined 
as the average size of finite clusters which includes a 
randomly-selected site, diverges at the transition. Their 
scaling properties, such as the exact critical exponents in 
two dimension Q and the upper critical dimensionality 
du = 6 [H, are well understood. 

Recently, the cluster heterogeneity H, defined as the 
number of distinct cluster sizes, is suggested as a useful 
indicator of a percolation transition . Consider a clus- 
ter size distribution function n s which is defined as the 
number of clusters of size s per site. When p = 0, the 
distribution is monodisperse with H = 1. As p increases, 
clusters nucleate and aggregate into larger ones. Conse- 
quently, in the subcritical phase, the cluster size distribu- 
tion becomes broader. In the supercritical phase, finite 
clusters are absorbed into the giant cluster to decrease H . 
Hence one expects that the cluster heterogeneity may be 
maximal at the transition. 

Lee et al. [|[ found that the maximum heterogeneity 
points indeed converge to a percolation threshold in the 
thermodynamic limit in the study of the so-called ex- 
plosive percolation Q. Exploiting the finite-size-scaling 
(FSS) property at the maximum heterogeneity points, 
they could clarify the nature of the explosive percolation 
transition. However, the scaling property of H by itself 
has not been understood fully even in the ordinary ran- 



dom percolations. In this paper, we establish the critical 
scaling law for H with emphasis on its FSS theory. 

We begin with a brief review for the scaling theory. 
For detailed review, we refer readers to Ref. p], Q . The 
ordinary percolation exhibits a continuous transition [f| . 
The correlation length diverges algebraically as £ ~ \z\~ v 
where e = p—p c and v is the correlation length exponent. 
At p = p c , the giant cluster is a fractal characterized with 
a fractal dimension df. Physical quantities have a singu- 
lar dependence on e. For example, the order parameter 
scales as m ~ for p > p c with the order parameter ex- 
ponent (3. The mean cluster size S diverges as S ~ |e|~ 7 
with the susceptibility exponent 7. The cluster size dis- 
tribution function n s {p) for finite s scales as 



n s (p) 



s/s 



(1) 



where t is called the Fisher exponent and s c is the char- 
acteristic cluster size. It diverges as s c ~ | e | 1 / with a 
critical exponent a. There exist scaling relations among 
those critical exponents. So, any exponent can be writ- 
ten in terms of two independent ones, say v and df. The 
scaling relations read /3 = v(d — df), 7 = v(2df — d), 
a = l/{vdf), and r = 1 + d/df. 

In finite systems of linear size L, the scaling laws are 
modified because the correlation length £ is limited by 
L. The FSS hypothesis assumes that a finite-size effect 
comes into play through the ratio between L and £ Q. 
It leads to the FSS ansatz for an observable Q(p, L) as 



Q{p,L) = L x QT Q {eL 1 '») , 



(2) 



where Xq is a scaling exponent and J : q{x) is a scaling 
function for Q Q. The FSS ansatz claims that a finite 
system suffers from a finite-size effect in the region |e| < 
i -1 /" whereas it behaves as an infinite one elsewhere. 
The aim of this paper is to find a scaling law for H. 

This paper is organized as follows: In Sec. [Hi we 
present numerical data for H in site and bond perco- 
lations in two-dimensional (2D) square and triangular 
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FIG. 1. (Color online) Cluster heterogeneity in site percola- 
tions in (a) and (c), and bond percolations in (b) and (d) on 
2D square lattices in (a) and (b) and triangular lattices in (c) 
and (d). The dotted lines represent the critical percolation 
threshold. Lattice sizes are L — 2 5 , ■ ■ • , 2 12 for the square lat- 
tices and L = 27 x 2° , ■ ■ ■ , 27 x 2 7 for the triangular lattices. 
The larger L is, the higher H is. 



lattices, which shows that the standard FSS form ([2]) is 
not valid for H . In Sec. IIIH we derive the central result 
of Eq. (j!5[) and present numerical data in d = 2, • • • , 6 
dimensional hypercubic lattices to confirm it. We sum- 
marize and conclude this paper in Sec. [IVJ 



II. PERCOLATIONS IN 2D 

We have performed Monte Carlo simulations of site 
and bond percolations in square and triangular lattices 
of L x L sites in 2D using the Newman and Ziff algo- 
rithm Q. This algorithm allows an efficient and fast 
measurement of a quantity Q as a function of the num- 
ber n of occupied sites or bonds. A corresponding quan- 
tity as a function of the occupation probability p is then 
obtained by the convolution Q(p) — ^ jn B{N,n 1 p)Q{n) 
where N is the total number of sites or bonds and 
B(N,n,p) = „, ( ^L), p"(l - p) N ~ n is the binomial dis- 
tribution function || Q . 

The cluster heterogeneity, averaged over N$ = 10 5 
samples, is presented in Fig. [TJ In all cases, the curves 
have a diverging peak (maximum cluster heterogeneity) 
at a position denoted by (p*,H*). The peak position 
seems to approach from below the critical point, which is 
p c = 0.592746 for the square lattice site percolation Q, 
1/2 for the square lattice bond percolation and the tri- 
angular lattice site percolation, and 2sin(-7r/18) for the 
triangular lattice bond percolation FSS properties 
are analyzed in Fig. [2j We find that (p c — p*) and H* 
scale algebraically with L with apparently universal ex- 
ponents. 

It is noteworthy that (p c — p*) does not scale as L^ 1 /" 
with the correlation length exponent v = 4/3 in 2D. Such 
a scaling would be natural from the FSS hypothesis of 
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FIG. 2. (Color online) FSS of (p c -p*) in (a) and H* in (b) for 
the site (open symbols) and the bond (filled symbols) percola- 
tions on 2D square (square symbols) and triangular (triangle 
symbols) lattices. All straight lines are guides to eyes. In (a), 
all the solid lines have the same slope while the dashed line 
has a different slope of 1/v = 3/4. 



Eq. ([2]). However, Fig. [2ja) excludes such a possibility 
definitely. It calls for an appropriate FSS theory for H . 



III. FSS THEORY FOR H 

In order to characterize the cluster heterogeneity, one 
needs to specify the size of each cluster. Jan et al. 
considered a FSS behavior of the rth largest cluster size 
at the critical point where the cluster size distribution 
follows a power law n s ~ s~ r . The rth cluster size s r is 
estimated from the relation 



r ~ L d I ds s" 7 



which yields that 

— 1 d 

s r ~ r T - 1 L t- 1 



(3) 



(4) 



We extend their idea to the off-critical region to derive 
the FSS theory for H. 

For p p ci the cluster size distribution function is 
given by Eq. ([T]) [l2|. Then, the size of a rth largest 
cluster is obtained from 



ds 



d„l-7 



L a s 



r(l - r, s r /s c ) , (5) 



where function T(u, x) = J^°dtt u 1 e ' is the incomplete 
gamma function. The characteristic size itself displays a 
FSS behavior 



| e | _ for | £ | > L -y v ^ 
L d f for lei < L~ x l v . 



(0) 



We now compare the average size of rth and (r + l)th 
clusters. Using Eq. ([S]), one finds that 



A r s = s r - s r+ i ~ L d s T e s / s " 



(7) 
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The distribution is dense when A r s < 1 and sparse when 
A r s > 1. The two regions are separated at s = Sq satis- 
fying 



L- d sle SQ l s " 



0(1) 



The rank of a cluster of size sq is given by 
r ~L d sl- T T(l-r,s /s c ) . 



(8) 



(9) 



Note that L d n SQ = 0(1). This implies that there are at 
least 0(1) clusters of all sizes s < sq. On the other hand, 
there are r clusters in the region s > Sq, whose sizes are 
distinct because A r s > 1. Therefore we find that 



H ~ s + r 



(10) 



To obtain the solution of Eq. ((SJ for so, we first assume 
that so *C s c . Then, the exponential term is negligible 
and the solution is given by 



so 



(11) 



Using the asymptotic behavior r(u, x — > 0) ~ — ^x u , one 
also finds that r ~ Sq. This solution is self-consistent 
when L d / T < kt 1/a or 



lei < L-V"fl 



with the exponent i/# 



r 



(12) 



(13) 



In the opposite case so ^ s c , the leading order solution 
of Eq. © is given by 



so 



s r In L 



(14) 



From the asymptotic behavior T(u,x — > oo) ~ e x , 
one also finds that tq ~ s c . This solution is self-consistent 
when |e| > L~ X I VH . 

Using the solution for so and tq , we can summarize the 
scaling property of H with the FSS form 

H(p, L) = L d l T T H {£L llvH ) • (15) 

The scaling function has a limiting behavior 

r d~ 1/<T In Id for \x\ > 1, 



constant for Ixl <S 1 



(16) 



so that Jf~ (\nL)\e\- 1 /' T for |e| » L" 1 /^ and if - i d / T 
for |e| < L~ X I VH . 

Remarkably, the FSS exponent vu for the cluster het- 
erogeneity is distinct from the correlation length expo- 
nent v. It explains the reason why (p c — p*) does not 
scales as L- X l y in Fig. [2ja). Instead, the numerical 
data are consistent with the scaling (p c — p*) ~ L~ X I VH 
with v H = 187/72 in 2D (see Table The solid lines 
in Fig. [Ha) have the slope of 1/vh- The maximum 
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FIG. 3. (Color online) Scaling analysis of H for the site ((a) 
and (c)) and bond ((b) and (d)) percolations on 2D square 
((a) and (b)) and triangular ((c) and (d)) lattices. The same 
data sets are used as in Fig. [1] 



heterogeneity is expected to scale as H* ~ L d l T with 
d/r = 182/187 in 2D. The solid line in Fig.^b) is of this 
slope and in agreement with the numerical data. The 
FSS scaling form is tested in Fig. [3l where we replot the 
data in Fig. [T] according to Eq. (|T5|) . All the data from 
different sizes collapse perfectly onto a single curve. 

We examine whether the FSS form in Eq. (|15j) for H 
is valid universally in higher dimensions up to the upper 
critical dimension d u = 6. We have performed exten- 
sive numerical simulations of site and bond percolations 
to measure H in hypercubic lattices in d = 2, • • • ,6 di- 
mensions. The system sizes are L = 2 3 ,---,2 8 in 3D, 
L = 2 3 , ■ ■ • , L 6 in 4D, L = 8, ■ ■ • , 28 with AL = 4 for 5D, 
and L = 6, • • • , 16 with AL = 2 for 6D. 

In Fig. [4ja) , (p c — p*) is plotted against 1/L for the 
site percolations. Also drawn are the straight lines corre- 
sponding to the scalings L 1 ^ and L 1 ^" . The parameter 
values used in the plot are collected from literatures and 
listed in Table HJ In all dimensions, (p c — p*) scales with 
the exponent 1/vh instead of 1/v. Figure E|b) shows 
the plot of H* against L, which also confirms the scaling 
H* ~ L d / T . We have obtained the same result for the 
bond percolations, whose data are not shown here. 



IV. SUMMARY AND DISCUSSION 

We have shown that the cluster heterogeneity diverges 
at the percolation critical point and that it satisfies the 
FSS form in Eq. |[ISJ). The FSS is governed not by the 
correlation length exponent v but by the new exponent 
i>h = (1 + df/d)v. Consequently, the maximum clus- 
ter heterogeneity points follow the scaling (p c —p*(L)) ~ 
j^-xjvH xhis is contrasted with the standard L~ x l v scal- 
ing of the effective critical points obtained from the order 
parameter data and the average cluster size data. 
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TABLE I. Critical points and critical exponents in the site percolations in d-dimensional hypercubic lattices. 

d = 2 d = 3 d = 4 d = 5 d = 6 

Pc (site) 0.592746 [8] 0.3116077 [13] 0.196889 [16] 0.140765 [17] 0.109017 [17] 

4/3 0.875 [14] 0.689 [15] 0.51 [18] 1/2 (MF) 

d/ 91/48 2.523 [14] 3.05 [16] 3.54 [16] 4 (MF) 

i/ir = (l + df/d)v 187/72 1.611 1.21 0.87 5/6 (MF) 

d/T = d/(l + d/d f ) 182/187 1.370 1.73 2.07 12/5 (MF) 
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FIG. 4. (Color online) (a) (p c — p*) vs. 1/L. The solid lines 
have a slope of 1/vh while the dotted lines have a slope of 
\jv. (b) H* vs. L with the solid lines having a slope of d/r. 



Since vu > v, the maximum heterogeneity points 
p*(L) constitute a singular path of e ~ —L~ X I VH in the 
subcritical phase to the critical point. Although they 
converge to p c in the L — > oo limit, the system along 
the path remains outside the critical region defined by 



the condition |e| <C L~ x l v . Physical quantities along a 
singular path may exhibit a peculiar FSS behavior. 

Consider, for example, the order parameter to which 
follows the FSS form in Eq. (0 with X m = -0/l>. In 
the subcritical phase, it scales with L as to ~ L~ d in the 
leading order [19[, which requires that the scaling func- 
tion should have a limiting behavior J- m (x <C —1) ~ 
\^ x \-(dv~p) _ g Q foe order parameter to*, evaluated at 
the maximum heterogeneity points, follows the scaling 
m* ~ L-< 3 l u J : m (-L- 1 l UH+1 / u ) ~ L-P/v-d s (i-v/v H )^ in _ 
stead of the standard critical scaling to* ~ L^ 13 ^. In 
fact, this is a manifestation of the crossover phenomena 
near a multicritical point in the context of the FSS [2(| ■ 
When there exist multiple relevant scaling fields for a 
critical point, scaling behaviors become dependent on a 
path approaching it, in general. Since 1/L is one of such 
relevant scaling variables 0,1111, one can naturally expect 
a nontrivial FSS along a singular path. The maximum 
cluster heterogeneity condition indeed allows us to access 
a singular path to the percolation critical point. 

This work was supported by Mid-career Researcher 
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